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1 Abstract 

A general proecological urban road traffic control idea for the street canyon is proposed with 
emphasis placed on development of advanced continuum field gasdynamical (hydrodynam- 
ical) control model of the street canyon. The continuum field model of optimal control of 
street canyon is studied. The mathematical physics approach (Eulerian approach) to vehic- 
ular movement, to pollutants' emission and to pollutants' dynamics is used. The rigorous 
mathematical model is presented, using gasdynamical (hydro dynamical) theory for both air 
constituents and vehicles, including many types of vehicles and many types of pollutant 
(exhaust gases) emitted from vehicles. The six optimal control problems are formulated and 
numerical simulations are performed. Comparison with measurements are provided. General 
traffic engineering conclusions are inferred. [PCCMM'99 poster session]. 

2 Description of the model 

The vehicular flow is multilane bidirectional one-level one-dimensional rectilinear and it is 
considered with two coordinated signalized junctions [1, 2]. The vehicles belong to different 
vehicular classes: passenger cars, trucks. The emission from vehicles are based on technical 
measurements and many types of pollutants are considered (carbon monoxide CO, hydro- 
carbons HC, nitrogen oxides NOx). The vehicular dynamics is based on hydrodynamical 
approach [3]. The governing equations are continuity equation for number of vehicles and 
Greenshields' equilibrium speed-density u-k model [4] . The model of dynamics of pollutants 
is also hydrodynamical. The model consists of a set of mutually interconnected vector non- 
linear, spatially three-dimensional, temporally dependent, partial differential equations with 
nonzero righthand sides (sources), and of boundary and initial problem. The pollutants, 
oxygen, and remaining gaseous constituents of air, are treated as mixture of noninteracting, 
Newtonian, viscous fluid (perfect or ideal gases). The model incorporates as variables the 
following fields: density of mixture, mass concentrations of constituents of mixture, velocity 
of mixture, temperature of mixture, pressure of mixture, and intrinsic (internal) energy of 
mixture. The model is based on assumption on local laws of balance (of conservation): of 
mass of mixture, of masses of constituents of mixture, of momentum of mixture, of energy 
of mixture, and of state equation (Clapeyron's law). The model of dynamics is solved by 
finite difference scheme. 
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The six optimization problems are formulated by defining functional of total travel 
time, of global emissions of pollutants, and of global concentrations of pollutants, both in 
the studied street canyon, and in its two nearest neighbour substitute canyons. The vector 
of control is a five-tuple composed of two cycle times, two green times, and one offset time 
between the traffic lights. The optimal control problem consists of minimization of the six 
functionals over admissible control domain. 



3 Equations of dynamics. 

El. Balance of momentum - Navier Stokes equation. 
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E2. Balance of mass - Equation of continuity. 



^ + div(pv) = S. (2) 



E3. Balance of mass of constituents - Diffusion equations. 
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E4. Balance of energy. 
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E5. Equation of state - Constitutive equation- Clapeyron's equation. 

p R 



p m, 



-T, (6) 



air 



Pi = Ci^-p,i = 1, ...,N. (7) 
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E6. Balance of vehicles - Equation of continuity. 
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+ div(^ t <„ t ) = 0. (9) 
vt = l,...,VT,l = l,...,n L ,r = l,...,n R . 

4 Optimization problems. 

E8. Control. 

u=(g 1 ,C 1 ,g 2 ,C 2 ,F)eU adm , (10) 

where <?j, are green times, Q, are cycle times, F is offset time, and U £Ldm is a set of admissible 
control variables. 

Fl. Total travel time for single canyon. 

J* TT = J TTT (u*) = M{u G f/ adm : Jttt(m)}, (11) 
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F2. Global emission for single canyon. 

J* = J E ( M *) = in f{ M G f/ adm : J E («)}, (13) 
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F3. Global pollutants' concentration for single canyon. 

J* c = J c (u*) = inf{u e £/ adm : J c (u)}, 



J c(u) = V/9stp / / / / Ci(x,y,z,i)dxdydzdt. 
~[ Jo Jo Jo Jo 

F4. Total travel time for canyon in street subnetwork. 
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F5. Global emission for canyon in street subnetwork. 

^E,ext = JeM^*) = mf{« e U adm : J E ,«t(«)}, 
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F6. Global pollutants' concentration for canyon in street subnetwork. 
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5 Conclusions. 



We assumed the data from real street canyon Krasihski Avenue in Cracow [5]. In general 
the numerical results are in very good agreement with measured data from [5]. 

In general some vehicular traffic and hydrodynamical parameters influence the solutions 
of optimization problems Fl, F2, F3, F4, F5, F6, but not all of them: 

Rl. The direction of velocity of air mixture is important. The optimal pollutant con- 
centrations for both single canyon and canyon with its two substitute nearest neighbour 
canyons, F3, F6, are the lowest if the velocity components of the boundary and initial value 
problems are equal to zero; further, they are the grater if velocity has only nonzero vertical 
component v z ; then, they are greater if velocity has only nonzero x-component v x ; further, 
they are again greater if velocity has two nonzero components (0, v y , v z ); next, they are again 
greater if velocity has three nonzero components (v x ,v y ,v z ); finally, they are the greatest if 
the velocity has only nonzero y-component v y . 

R2. The optimal values for Fl, F4 cases, for F2, F5 cases, and for F3, F6 cases, 
decrease with "uniformization" of vehicles, when we pass from nonuniform vehicles to uniform 
ones. It is a result of decrement of the number of vehicles moving in the canyon. For uniform 
the values of maximum free flow speed, jam, saturation, threshold, green, and red densities 
take on minima. 

R3. The long vehicular queues decrease total travel times Fl, F4, and they increase 
both optimal emissions F2, F5, and concentrations of pollutants F3, F6. The decrement 
of total travel times Fl, F4, with long vehicular queues is result of clustering of vehicles. 

R4. The constant of temperature scale does not differentiate the values of optimal 
concentrations of pollutants F3, F6 in the temperature range near standard temperature 
and pressure STP conditions. However, it diminishes them even hundredfold for very high 
temperatures. 

R5. The functional form of initial and boundary conditions affects the optima. If they 
are constant then optimal concentrations F3, F6 are twice higher than in the case when 
they are changing exponentially in space in three dimensions. 

R6. The presence of vehicles on both left and right lanes is important. The optimal 
total travel times and emissions are halved in absence of vehicles on left or right lanes with 
respect to situation when they circulate on both left and right lanes. 

R7. The values of saturation, arrival, or jam vehicular density, and of vehicular free flow 
velocities also affect the optima Fl, F2, F3, F4, F5, F6. 

R8. The assumption of energy conservation equation, of t her mo diffusion effect, of chem- 
ical potential and of Grand Canonical ensemble, as well as of influence of gravity on intrinsic 
energy and on chemical potential, drastically changes the optimal concentrations F3, F6 
towards measured ones [2]. 

R9. The value of time of simulation and of discretization in time affects much the 
optima. The values of optimal solutions Fl, F2, F3, F4, F5, F6 increase from tenfold to 
hundredfold. Also the optimal 5-tuples for Fl, F2, F3, F4, F5, F6 change their values. It 
is result of cumulative effect of length of period of simulation T on integral functionals Fl, 
F2, F3, F4, F5, F6. 
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